GLOBAL SOLUTION AND TIME DECAY OF THE COMPRESSIBLE 
EULER-MAXWELL SYSTEM IN R 3 
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Abstract. We study the global existence and large time behavior of solutions near a constant 
equilibrium to the compressible Euler-Maxwell system in R 3 . The previous works mainly as- 
sumed that the H (N > 3) and L 1 norms of the initial data are sufficiently small. In this 
paper, we first construct the global unique solution by assuming that the H 3 norm of the initial 
data is small, but the higher order derivatives can be arbitrarily large. If further the initial data 
belongs to H~ s (0 < s < 3/2) or 13^ oo (0 < s < 3/2), by a regularity interpolation trick, we 
obtain the various decay rates of the solution and its higher order derivatives. As an immediate 
byproduct, the L p -L 2 (1 < p < 2) type of the decay rates follow without requiring that the L p 
norm of initial data is small. 



1. Introduction 

The dynamics of the electrons interacting with their self-consistent electromagnetic field can 
be described by the compressible Euler-Maxwell system [12] : 

' dth + div(nu) = 0, 

dt{nu) + div(nn (£> u) + Vp(n) = —fi{E + eu x B) — -nu, 

eX 2 d t E - Vx5 = ehu, , . 

ed t B + V x E = 0, ^ ' ' 

\ 2 &i\E = TT-oo — h, divB = 0, 
. (h,u,E,B)\ t=0 = (h ,u ,E ,B ). 

The unknown functions n, u, E, B represent the electron density, electron velocity, electric field 
and magnetic field, respectively. We assume the pressure p(h) = An 1 with constants A > and 
7 > 1 the adiabatic exponent. 1/r > is the relaxation time. A > is the Debye length, and 
e = 1/c with c the speed of light. In the motion of the fluid, due to the greater inertia the ions 
merely provide a constant charged background > 0. 

Despite its physical importance, due to the complexity there are only few mathematical 
studies on the Euler-Maxwell system. In one space dimension, Chen, Jerome and Wang [1, 
proved the global existence of entropy weak solutions to the initial-boundary value problem 
for arbitrarily large initial data in L°° . Since the Euler-Maxwell system is a symmetrizable 
hyperbolic system, the Cauchy problem in M. 3 has a local unique smooth solution when the 
initial data is smooth, see Kato |1U] and Jerome [U] for instance. Recently, there are some 
results on the global existence and the large time behavior of smooth solutions with small 
perturbations, see Duan [2], Ueda and Kawashima Ueda, Wang and Kawashima [22J. For 
the asymptotic limits that derive simplified models starting from the Euler-Maxwell system, we 
refer to [5J [T51 12"5] for the relaxation limit; for the non-relativistic limit; [1'6\ [T5] for the 
quasi-neutral limit; [19} [20] for WKB asymptotics; and references therein. 

The main purposes of this paper are to refine a global existence of smooth solutions near the 
constant equilibrium (rioo, 0, 0, i?oo) to the compressible Euler-Maxwell system and to derive 
some various time decay rates of the solution as well as its spatial derivatives of any order. 

2010 Mathematics Subject Classification. 83C22; 82D37; 76N10; 35Q35; 35B40. 

Key words and phrases. Compressible Euler-Maxwell system; Global solution; Time decay rate; Energy 
method; Interpolation. 

Corresponding author: Yong Wang, wangyongxmu@163.com. 

Supported by the National Natural Science Foundation of China-NSAF (No. 10976026) and the Natural 
Science Foundation of Fujian Province of China (No. 2012J05011). 



2 



ZHONG TAN, YANJIN WANG, AND YONG WANG 



We should emphasize that our results highly rely on that we consider the relaxation case. The 
non-relaxation case is much more difficult, we refer to [31 [5] for such direction. It turns out that 
it is more convenient to reformulate the compressible Euler-Maxwell system (jl.ip as follows. 
Without loss of generality, we take the constants r, s, A, A, in (jl.ip to be one. We define 



n(x, t) 



2 

7-1 



7-1 
2 



1 



u(x, t) 



E(x,t) = ±E(x,^, B(x,t) = ^B(a 

Then the Euler-Maxwell system (jl.ip is reformulated equivalently as 

' dtn + divu = — u ■ Vn — /undivu, 
dtu + vu + u x Boo + Vn + vE = —u ■ Vit — jinSJn — u x B, 
dtE — uV x B — uu = uf(n)u, 
d t B + z/V x E = 0, 
div£ = -vf(n), divB = 0, 
(n,u,E,B)\ t=0 = (no,u ,E ,B ). 

:= and the nonlinear function f(n) is defined by 

2 

1 



(1.2) 



(1.3) 



Here [i := ^n— , 



1 + 



7' 



-n 



1. 



(1.4) 



(1.5) 



Notice that we have assumed 7 > 1. If 7 = 1, we instead define 

n := y/A (In h — In Uqo) = V^Aln n. 

In this paper, we only consider the case 7 > 1, and the case 7 = 1 can be treated in the same 
way by using the reformulation in terms of the new variables correspondingly. 

Notation. In this paper, we use H S (R 3 ), s G M to denote the usual Sobolev spaces with norm 
\\-\\jjs and L P (M 3 ),1 < p < 00 to denote the usual LP spaces with norm ||-|| LP . V e with an 
integer £ > stands for the usual any spatial derivatives of order t. When £ < or £ is not a 
positive integer, stands for defined by A / := ^" _1 (|f |^^/), where & is the usual Fourier 
transform operator and is its inverse. We use H s (U. 3 ),s £ f to denote the homogeneous 
Sobolev spaces on M 3 with norm \\-\\jj 3 defined by H/H^s := ||A S /|| L2 . We then recall the 
homogeneous Besov spaces. Let (j) £ C£°(J?f) be such that 0(f) = 1 when |f| < 1 and 0(f) = 
when |f | > 2. Let 93(f) = 0(f ) — 0(2f ) and <fj(£ ) = (/'(2~' J f ) for j € Z. Then by the construction, 
Yl PjiO = 1 if f / 0. We define Ajf := ^~ 1 ((fj)*f, then for s G M. we define the homogeneous 



Besov spaces i?f ^(M 3 ) with norm 



B:) 



defined by 



sup2^ 



L- 



(1.6) 



Throughout this paper we let C denote some positive (generally large) universal constants 
and A denote some positive (generally small) universal constants. They do not depend on either 
k or N; otherwise, we will denote them by Cjv, etc. We will use a < b if a < Cb, and a ~ b 
means that a < b and b < a. We use Co to denote the constants depending on the initial data 
and k,N,s. For simplicity, we write -B)||x := ll-^llx + H-^H x an d I f := Im 3 f dx. 

For N > 3, we define the energy functional by 

N 



z=o 



2 

L 2 



(1.7) 



and the corresponding dissipation rate by 

N 

V N (t) :=^2\\v l (n,u) 



1=0 



N-l 

1=0 



WE 



AT-1 

i=i 



V L B 



2 

L 2 



(I- 
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Our first main result about the global unique solution to the system (jl.3p is stated as follows. 
Theorem 1.1. Assume the initial data satisfy the compatible conditions 

div£ = -vf{no), divB = 0. (1.9) 
There exists a sufficiently small 5q > such that if £$(0) < Sq, then there exists a unique global 



solution (n,u, E, B)(t) to the Euler- Maxwell system (|1.3|) satisfying 

roc 

sup £ 3 (t)+ / V 3 (r)dT < C£ 3 (0). 

0<t<oo JO 



(1.10) 



Furthermore, if £n(0) < +oo for any N > 3, there exists an increasing continuous function 
Pn(') with P/v(0) = such that the unique solution satisfies 



poo 

sup £N{t) + j V n (t) dr < P N (£ N 
0<t<oo Jo 



(0)). 



(1.11) 



The proof of Theorem 11.11 is inspired by the recent works of Guo [6] and Wang [23J . The 
major difficulty here is the regularity-loss of the electromagnetic field. We will do the refined 
energy estimates stated in Lemma which allow us to deduce 

d :£s + V 3 < (1-12) 



and for N > 4, 



d_ 

dt 



dt 



£n + T^n < CnT>n-\£n ■ 



(1.13) 



Then Theorem 11.11 follows in the fashion of [61 [23] . 

Our second main result is on some various decay rates of the solution to the system (|1 .3j) by 
making the much stronger assumption on the initial data. 

Theorem 1.2. Assume that (n,u, E, B)(t) is the solution to the Euler-Maxwell system (jl.3p 
constructed in Theorem \l.l\ with N > 5. There exists a sufficiently small 5o = o~o(N) such that if 
£n{0) < So, and assuming that (uo, Eq, Bq) G H~ s for some s £ [0,3/2) or (uo, Eq, Bq) G B^^ 
for some s G (0,3/2] ; then we have 

\\(u,E,B)<t)\\ A -.<C (1.14) 

or 

\\{u,E,B){t)\\ K s <Co. (1.15) 
Moreover, for any fixed integer k > 0, if N > 2k + 2 + s, then 

<Co(l + t)~ — . 



V k (n,u,E,B)(t) 



I? 



Furthermore, for any fixed integer k > 0, if N > 2k + 4 + s, then 

V k {n,u,E){t) 



L 2 



<C7 (l + r^: 



if N > 2k + 6 + s, then 
if N > 2k + 10 + s and B a 



V k n{t) 
; then 



L 1 



k + 2+s 

<C (l + t) — : 



V k (n,dwu)(t) <C (l + ty^ + i +s \ 



L 2 



(1.16) 

(1.17) 
(1.18) 
(1.19) 



The proof of Theorem 1 1.2 1 is based on the regularity interpolation method developed in Strain 
and Guo [T7], Guo and Wang [7] and Sohinger and Strain [16]. To prove the optimal decay 
rate of the dissipative equations in the whole space, Guo and Wang [?] developed a general 
energy method of using a family of scaled energy estimates with minimum derivative counts 
and interpolations among them. Note that the homogeneous Sobolev space H~ s was introduced 
there to enhance the decay rates. By the usual embedding theorem, we know that for p G (1, 2], 
LP C H~ s with s = 3(| - |) G [0,3/2). Hence the LP-L? type of the optimal decay results 
follows as a corollary. However, this does not cover the case p = 1. To amend this, Sohinger 
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and Strain |16j instead introduced the homogeneous Besov space B 2 ^ due to the fact that the 



endpoint embedding L C B 2 ^ holds. The method of [3[T6] can be applied to many dissipative 
equations in the whole space, however, it can not be applied directly to the compressible Euler- 
Maxwell system which is of regularity-loss. Based on the refined energy estimates stated in 
Lemma [2THT42.91 we deduce 



dt k 
lfc+2 



+ v k+2 <c k 



n, u 



V 



k+2. 



n,u 



L 2 



V k+2 (E,B) 



L 2 



(1.20) 



where £% +2 and V k k +2 with minimum derivative counts are defined by (|3.9p and (|3.10|) respec- 
tively. Then combining the methods of [16] and a trick of Strain and Guo [T7] to treat the 
electromagnetic field, we are able to conclude the decay rate f)1.16j) . If in view of the whole solu- 
tion, the decay rate (|1.16p can be regarded as be optimal. The faster decay rates ()1.17p - (|1.19p 
follow by revisiting the equations carefully. In particular, we will use a bootstrap argument to 
derive (fl~T9|) . 

As quoted above, by Theorem II. 2\ we have the following corollary of the usual L p -L 2 type 
of the decay results: 

Corollary 1.3. Under the assumptions of Theorem \1.2\ except that we replace the H~ s or B^^ 

assumption by that (uo, Eq, Bq) £ L p for some p G [1,2], then for any fixed integer k > 0, if 
N > 2k + 2 + s p , then 



Here the number s 



V k (n,u,E,B)(t) <C (l + t) 



k + 3p 
2 



L 2 



(1.21) 



v ■" 



MM 



Furthermore, for any fixed integer k > 0, if N > 2k + 4 + s p , then 



V k (n,u,E)(t) <C Q (l + i) 



v 

k + l + s p 



L 2 



ifN>2k + 6 + s p , then 



if N > 2k + 10 + s p and B 



V k n(t) 
0, then 



L 2 



<C (l + t) r 



V k (n,divu)(t) <C (l + t)-^ + l 



+ Sp) 



L 2 



1.221 



(1.23) 



(1.24) 



The followings are several remarks for Theorem ll.il Theorem 1 1 . 2 1 and Corollary 11.31 

Remark 1.4. In Theorem M.R we only assume the H 3 norm of the initial data is small, but the 
higher order derivatives can be arbitrarily large. Notice that in Theorem the H~ s and B^'^ 
norms of the solution are preserved along the time evolution, however, in Corollary \1.3\ it is 
difficult to show that the LP norm of the solution can be preserved. Note that the L? decay rate 
of the higher order spatial derivatives of the solution are obtained. Then the general optimal L q 
(2 < q < oo) decay rates of the solution follow by the Sobolev interpolation. 

Remark 1.5. We remark that Corollary \Tj^ not only provides an alternative approach to derive 
the LP-L 2 type of the optimal decay results but also improves the previous results of the LP-L 2 
approach in Ueda and Kawashima [21] and Duan [2]. Ln Ueda and Kawashima |21j . the decay 
rates (jl.2ip - (jl.23p with p = 2 were proved by using the time weighted energy method, and when 
p = 1 they were proved by combining the time weighted energy method and the linear decay 
analysis but under the stronger assumption that \\(uq, uq, Eq, Bq)\\ l i is sufficiently small. Ln 
Duan [2\, assuming that B^ = and \\(uq, Eq, Bq)\\ l i is sufficiently small, by combining the 
energy method and the linear decay analysis, Duan proved that 

\\n(t)\\ L2 <C (l + t)-^, \\(u,E)(t)\\ L2 <C (l + t)-i and \\B(t)\\ L2 < C (l + t)~i. (1.25) 

Notice that for p = 1, our decay rate of n(t) in (|1.24p is (1 + t) -13//4 . 

The rest of our paper is organized as follows. In section [21 we establish the refined energy 
estimates for the solution and derive the negative Sobolev and Besov estimates. Theorem 11.11 
and Theorem 11.21 are proved in section [3[ 
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2. Nonlinear energy estimates 



In this section, we will do the a priori estimate by assuming that H 3 < i <C 1. Recall 

the expression (jl.4p of /(n). Then by Taylor's formula and Sobolev's inequality, we have 



/(n) ~ n and / (fe) (n) 



< Cfc for any k > 1. 



(2.1) 



2.1. Preliminary. In this subsection, we collect the analytic tools which will be used in the 
paper and prove a basic estimate for the nonlinear function f(n). 



L 2 



Lemma 2.1. Let 2 < p < +oo and a,m,£> 0. Then we have 

||v a /|| LP <c p ||v m /||^||vV 

Here < 9 < 1 (if p = +oo, i/ien we require that < 9 < 1) and a satisfy 

Proof. For the case 2 < p < +oo, we refer to Lemma A.l in [7J; for the case p = +oo, we refer 
to Exercise 6.1.2 in g| (pp. 421). □ 



(2.2) 
(2.3) 



Lemma 2.2. For any integer k > 0, we have 



V k f(n) 



and 



<C k 

W L 2 



V fc n 



3/4 
L 2 



<c fc 



V k n 



L 2 



(2.4) 
(2.5) 



Proof. The proof is based on Lemma 12.11 For (|2.4p . we refer to Lemma 3.1 in [JJ. For (|2.5p . 
in light of ()2.ip . it suffices to prove that when fe > 1, (|2.5p holds for all f(n) with bounded 
derivatives. We will use an induction on k > 1. If /c = 1, we have 

l|V/(n)|| L2 H|/'^)Vn|| L2 <l|Vn|| L2 . 

Assume ()2.5p holds for from 1 to k — 1. We use the Leibniz formula to have 

V fc /(n) j2 

L 2 w 

/ k-1 



(2.6) 



V fc - 1 (/'(«)Vn 
/'(n)V fc n 



L 2 



+ 



V/ / (n)V fc - 1 n + ^ ||vV(n)V* 



71 



(2.7) 



L 2 



Here if k = 2, then the summing term in (j2.7j) is nothing, etc. By Holder's inequality and 
Sobolev's inequality, we have 



f'{n)V k n 



L 2 



+ 



Vf'{n)V k - l n 



< 



L 2 



V k n 



L 2 



+ ||Vn|| L3 





< 


V k n 






L 6 ~ 




L 2 



(21 



For the summing term we use the induction hypothesis to obtain that for 2 < £ < k — 1, 



VV(n)V fc ^n 
By Lemma |2~7L| if ^ < [^p] , then we have 



n 


< 


VV'(n) 




V k - e n 


< 


V e n 




V k ~ e n 




L 2 ~ 


L 2 








L 2 









V k ~ e n 




L 2 





< liv c 



n 



\l 2 













V k n 


\\n\\ T 2 
L 2 11 "L 2 


V k n 


* < 

L 2 


V k n 



L 2 



where a is defined by 



k-i+l f 
£ = a x ; + k x 



3A- 



A- 



a 



2(fc-£) + 3 



<3; 



if i> [*fl] + 1, then we have 







V k - £ n 




L 2 





Z \H L 2 



V k n 



|V°n||| 2 



V k n 



L 2 



< 



V k n 



L 2 



(2.9) 

(2.10) 
(2.11) 

(2.12) 
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where a is defined by 



k — £+-=ax-+kx 1 — — 
2 k \ k 



3k 

a = Ye <3 - 



We thus conclude the lemma. 

We recall the following commutator estimate: 
Lemma 2.3. Let k > 1 be an integer and define the commutator 

V k ,g]h = V k (gh)-gV k h. 

Then we have 



(2.13) 
□ 



V fc ,<? 



h 



L 2 



<C k [\\Vg\\ La 





+ 


V k g 


J 






L 2 







(2.14) 
(2.15) 



Proof. It can be proved by using Lemma 12. 11 see Lemma 3.4 in [TT] (pp. 98) for instance. □ 

We have the LP embeddings: 
Lemma 2.4. Let < s < 3/2, 1 < p < 2 with 1/2 + s/3 = f/ien 

11/11*-. SII/IU,- (2.16) 

Proof. It follows from the Hardy-Littlewood-Sobolev theorem, see [3]. □ 
Lemma 2.5. £ei < s < 3/2, 1 < p < 2 wit/i 1/2 + s/3 = l/p, then 



B~ s ~ II-' Hi'' 

2 , OG 



Proof. See Lemma 4.6 in [16J. 

It is important to use the following special interpolation estimates: 
Lemma 2.6. Let s > and £ > 0, i/ien we /iai>e 



(2.17) 
□ 



< 



L 2 



l-6> 



L 2 



where 9 



1 



+ l + s 



Proof. It follows directly by the Parseval theorem and Holder's inequality. 
Lemma 2.7. Zei s > and £ > 0, i/ien we Ziaue 

1 



V7 



< 



L 2 



L 2 



^2,00 



where 9 



e + i + s 



Proof. See Lemma 4.5 in |16j . 



(2.18) 
□ 

(2.19) 
□ 



2.2. Energy estimates. In this subsection, we will derive the basic energy estimates for the 
solution to the Euler-Maxwell system (|1 .3j) . We begin with the standard energy estimates. 

Lemma 2.8. For any integer k > 0, we have 

r k+2 „ k+2 

d 
dt 



l=k 

< C k 1 1 hi. m 
+ \\(n,u) 



Y,\\v l (n,u,E,B) ' 2 + A^||v 



l=k 



2 

L 2 



'k+2 



+ \\VB\\ L2 ) (£||v'(n,u) L2 + E|| V ^ 



V 



k+2, 



n.u, 



L 2 



V k+2 (E,B 



l=k 



fc+i 



l=k 



L 2 



+ 



V k+1 B 



L 2 



(2.20) 
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Proof. The standard V' (I = k, k + 1, k + 2) energy estimates on the system (jl.3p yield 



1 d 

Ydt 



V(n,u,E,B) 



+ v 



Vu 



2 

L 2 



= -^J V(ndivu)Vn + V'(nVn) • V'u - ^ V(u ■ Vn)Vn + V{u- Vu) ■ Vu 

- J V l (u x B) ■ Vu + v J V(f(n)u) ■ VE 
:=h + h + h + h- 

We now estimate I\ ~ Z4. First, we use the commutator notation (|2.14p to rewrite I\ as 



V\n 



divu Vn + nV Vn + 



V\n 



divuVn + 



V',n 



Jl = — yu y (ndivVu + 

= -fi J ndiv(VuVn) 
By integrating by parts, we have 

- y ndiv(VuVn) = J Vn- VuVn < \\Vn\\ L 
We use the commutator estimate of Lemma 12.31 to bound 



V,n 



Vn ) • V l u 



Vn • Vn. 



(2.22) 



Vu 




Vn 






L 2 




L 2 



V,n 



divuVn < Q ( ||Vn|| io 



V 1 divu 


+ 


Vn 






L 2 




L 2 



||divu||j 



Vn 



L 2 



< Q ||V(n,u)|| ioo V'(n,«) 



L 2 



Vn 



(2.23) 



(2.24) 



L 2 



Bounding the last term of I\ similarly, and then applying the same arguments to the term I2, 
by Sobolev's and Cauchy's inequalities, we deduce 

2 



/l + /2<C,||(n,u 



V{n,u) 



L 2 



(2.25) 



Next, we estimate the term ^3, and we must be much more careful with this term since the 
magnetic field B has the weakest dissipative estimates. First of all, we have 

l 1 

(2.26) 

1=1 

Here, if I < 1, then it's nothing, and etc. We have to distinct the arguments by the value of I. 
First, let I = k. We take L 3 — L 6 and then apply Lemma 12. II to have 



h = -J2ctj V~ e u x Vb ■ Vu <c t J2 II V^uVb 



V k ~ e uV e B 


< 


V k - e u 




Vb 






L 2 ~ 




I? 







I 

II • U \\ L 2 



Vu 



L 2 



l|VS||^ 2 fc 



k 

L 2 



where a is defined by 



, „ 1 I , ( , * 

k — £+-=ax-+kx 1 — - 

2 k \ } 

Hence by Young's inequality, we have that for I = k, 



k k 



h<C k [ \\u\\ H k + || VB|| i2 



Vu 



L 2 



+ 



V +1 B 



(2.27) 



(2.28) 



(2.29) 



We then let I = k + 1. If 1 < I < A;, we take L 3 — L 6 and by Lemma 12 . 1 1 again to obtain 

(2.30) 



v k+l-i uv i B 


< 






Vs 






L 2 ~ 




L3 







7 Q 7;II fc „ 
~~ 11 • "11^,2 



V fc+1 u 



L 2 



1 '|VB||^ T 
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where a is defined by 



1 



k + l- £ + - = ax- + (k + l)x 1 



2 k 
if £ = k + 1, we take L°° - L 2 to get 



uV k+1 B 



k k 



< 



L- 



L- 



We thus have that for I = k + 1, by Sobolev's inequality, 



^3 < Cfc ( Hull fc 



+ HV5|| L2 



> 




2 

+ 


V k+1 B 








L 2 







We now let I = k + 2. If 1 < I < k, we take L 3 - L 6 and by LemmaO again to have 



v k+2-e uV e B 


< 


V k+2 ~ l u 










L 2 ~ 




L3 




L6 



||V Q u| 



V fc+2 u 



i-i 



1-4 



l|VB||^ 



L 2 



where a is defined by 

1 I 
- = a x - 
2 

if ^ = A; + 1 or k + 2, we take L°° - L 2 to get 



k + 2-£ + - = ax- + (k + 2)x[l--\^a = 2 + —<- + 2; 



k k 



VuV k+1 B 
and 

uV k+2 B 
We thus have that for I = k + 2, 



< 



L 2 



L- 



< 



V fc+2 fi 



L 2 



Is < Q (IMI fl4+w + [IVBII^ 
+ C||«|| £0 



> 


V fc+ \ 


2 

+ 


V fe+1 B 






L 2 









V fc+2 u 






L 2 







(2.31) 
(2.32) 
(2.33) 

(2.34) 

(2.35) 

(2.36) 
(2.37) 

(2.38) 



We now estimate the last term I4. We again have to distinct the arguments by the value of 
I. First, for I = k or k + 1, we have 

I 1 

(2.39) 

1 ; - 1 - 

If < £ < I - 2, we take L°° - L 2 and by Lemma |2~T1 and the estimate (|2.4f) of Lemma [2721 to 
obtain 

V £ /(")V^ £ n 



I r 1 

^=0 £=0 



< 



L 2 



V7(n) 



L 2 





1 




3 






4 


V f+2 n 


4 






L 2 




L 2 





<C; ( IMIl 2 ' 



1 21+3 



L ' 2 , 

21+3 
21 

L 2 



n\ 



L- 

II 2 

21+3 





|±2\ 4 















(2.40) 



L 2 



iiv a «n x r 



21 +:', 



L 2 



where a is defined by 



l-f=ax 



21 + 3 
21 



+ 1 x 1 



21 + 3 
2/ 



3/ 



a 



2£ + 3 



</; 



if £ = Z — 1 , we take L 6 — L 3 and by the estimate (|2.5p of Lemma 12.21 to have 



V ,-1 /(n)V« 



< 



L 2 



V'-7(«) L6 l|Vn|| x3 < Q 



L 2 



\H 2 ! 



(2.41) 
(2.42) 
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if I = I, we take L 2 — L°° and by the estimate (|2.4p of Lemma 12.21 to have 



L 2 



< 



V'/(n)« 

We thus have that for i = k or + 1, 

h < Q \\(n,u)\\ H i nH 2 
Now for I = k + 2, we rewrite as 

k+2 



L 2 





L 2 ' 







H 2 ■ 



V l (n,u) 



L 2 



+ 



V l E 



L 2 



(2.43) 



(2.44) 



k+2-i a , vjk+2 



f(n)V k+2 u + V k+2 f(n)u) • V fc+2 £ 



u ■ V K+Z E 



yk+2 j 



+ £ ^+2 / V (v fc + 2 -V(n)V^u) • V fc+1 £ 



f(n)V k+2 u + V k+2 f(n)u) -V K+ ' Z E 



ik+2 j 



(2.45) 



+ (A; + 2) y (v w /(n)VH + V/(n)V t+2 H ).V i+1 fi 

+ X5 C 'fc+2 / V k+3 - e f(n)V e u-V k+1 E + ^Ci +2 f V k+2 - e f(n)V e+1 u-V k+1 E 

O " D 1 " 



V k+2 u 



L 2 



+ 



V k+2 f(n] 



:— hi + I 42 + ^43- 
By Lemma 12.21 we have 

/4i<C fc (||/(n)|| z 

<C k \\(n,u)\\ Lc 

and 

h2 < C k (||v fc+2 /(n)|| L2 \\Vu\\ LOB + ||V/(n)|| 
<C k \\V(n,u)\\ L 



L 2 



\U\\ 



V k+2 E 



L 2 



V k+2 (n,u) 




yk+2 E 






L 2 




L 2 



(2.46) 



V k+2 u 














L 2 



V fc+2 (n,u) 










L 2 




L 2 



(2.47) 



As for the cases / = k, k + 1 for I4, we can bound ^43 by 



^43 < C fc ||u||^fc+i nH 2 ( V fe+1 (n,u) 
Hence, we have that for I = k + 2 



L 2 



+ 



V fc+1 £ 



74 < Ci, Ik 



V fc+1 (n,u) 



+ C fc || (n,u) || L oo V fc+2 (n,u) 



L 2 



L 2 



+ 



V k+2 (n,u) 



L 2 



+ 



2 

L 2 
V k+1 E 



(2.48) 



L 2 



V k+2 E 



(2.49) 



L 2 



Consequently, plugging the estimates for I\ ~ I 4 into (|2.2ip with I = k, k + 1, k + 2, and then 
summing up, we deduce (|2.2U|) . □ 

Note that in Lemma 12.81 we only derive the dissipative estimate of u. We now recover the 
dissipative estimates of n, E and B by constructing some interactive energy functionals in the 
following lemma. 
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Lemma 2.9. For any integer k > 0, we have that for any small fixed r/ > ; 



ft E / WwW]T Vu-V'E-n 
\i=k ^ i=k^ ^ 



V k E ■ V fc V x B 



/k+2 




k+i 


2 






E 

\l=k 


\w l n 


+y \\ve 

L 2 ^— ' II 

l=k 


+ 

L 2 




1.) 



fc+2 



Z=fc 



^ C E || V ' n [ 2 + C * ( II ( n ' u ) II + [| VB|& ) ( £ || V z (n, U ) 



'fc+2 



d=k 



L 2 



+ 



V k+1 B 



(2.50) 



Proof. We divide the proof into several steps. 
1: Dissipative estimate ofn. 
Applying V' (I = k, k + 1) to ()1.3[) 2 and then taking the L 2 inner product with VV'n, we 
obtain 



d t V l u ■ VV l n + 





2 




VV ; n 


< -1/ 




L 2 






■ 


V 



V'u 




V ?+1 n 






L 2 




L 2 



V* (u • Vn + ,unVn + n x 1?) 



L 2 



V m n 



(2.51) 



L 2 



The delicate first term on the left-hand side of (|2.5ip involves 9jV'n, and the key idea is to 
integrate by parts in the t- variable and use the continuity equation (|1 .3j) . Thus integrating by 
parts for both the t- and x- variables, we obtain 



J V l d t u ■ VV'n = f J V l u ■ VV're -J Vu- VV l d t n 



d 
d 



Vn • W'n 



V'divn 



> — / V l n • VV'ra - C 
~ dt 



V m n 



L 2 
2 



J V'divnV' (n ■ Vn + /zndivn) 



L 2 



c 



V(u ■ Vn) 



L 2 



c 



V z (ndivu) 



L 2 

(2.52) 



Using the commutator estimate of Lemma 12.31 we have 



V l (u ■ Vn) 



< 



L 2 



u • V'Vn 



< u 



L 2 

V l+1 n 



+ 



■ Vn 



L 2 



L 2 



+ Ci \\Vu\ 



V'n 



L 2 



+ c. 



L 2 



|Vn|| 



(2.53) 



<a||(n,«)||^3 



( 


V l (n,u) 


+ 


V m n 








L 2 







Similarly, 



Hence, we obtain 



V / (ndivn)| L2 <C / ||(n,n)|| / ,3( 



V\n,u) 


+ 


V l+1 u 






L 2 







(2.54) 



d 



dt 



V l d t u ■ VV'ra > -f / V l u ■ V'Vn - C 



Ci\\(n,u)\\ 



H"' 



V l (n,u) 



L 2 



2 

L 2 
+ 



V z+1 (n,n) 



(2.55) 
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Next, integrating by parts and using the equation (jl.3p 5 . we have 
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- J V l E ■ VV'n = v J V l divEV l n = -v 2 J V l f(n)V l n 



v 2 / V l (n + f(n) -n)V l n 



(2.56) 



Vn 



L 2 



+ Ci \\n\\ H 3 



2 

L 2 



Here we have used the estimate || V'(/(n) — ^)|| L 2 < Ci IMIh 3 II^HL 2 ' which follows by notic- 
ing that f(n) — n ~ n 2 and the similar arguments presented in Lemma 12.21 
Lastly, as in (|2.53p ~ (j2.54p . we have 



V (u ■ Vu + jLtnVn) 



L 2 



<C,||(n,u) 



\H 3 



V l (n,u) 


+ 


X7 l+1 (n,u) 






L 2 







From the estimate of 1^ in Lemma 12.81 we have that for I = k or k + 1 



V* (u x B) 



L 2 



<C k (\\u 



'H? +1 nH 2 



+ l|V5|| L2 



)( 




+ 


V k+1 B 


,) 






L 2 







(2.57) 
(2.58) 



d_ 



Plugging the estimates (|2.55[) — (|2.58[) into ()2.5ip . by Cauchy's inequality, we obtain 

fe+i „ k+2 



2 
L 2 



[ Vn- VV'n + A^ljv^ 

l=k J l=k 

|| V H[ 2 + ^ (ll^^ll^xnfl. + U^fe) ( E ||v 



n.u. 



L 2 



+ 



(2.59) 



This completes the dissipative estimate for n. 
Step 2: Dissipative estimate of E. 

Applying V 1 (I = k,k + 1) to (jl.3P a and then taking the L 2 inner product with V l E, we 
obtain 

2 



V'flfit • V l E + v 



X7 E 



L 2 



< -v J W l n-V l E + C 



Vu 




V l E 






L 2 




L 2 



+ 



V* (u • Vn + fmVn + u x B) 



L 2 



V l E 



(2.60) 



L 2 



Again, the delicate first term on the left-hand side of (|2.60p involves u, and the key idea 
is to integrate by parts in the t- variable and use the equation (|1.3P n in the Maxwell system. 
Thus we obtain 



(2.61) 



~ dt 

From the estimates of I4 in Lemma |2.8| we have that 



VE-J 


V l u 


V l d t E 








VE-v 


Vu 


W 









v J VV V* (/(n)u + V x B). 



V (f(n)u) 



L 2 



< C, \\u 



1 \\ a \\H l r\H 2 



L 2 



(2.62) 



We must be much more careful with the remaining term in (|2,6ip since there is no small factor 
in front of it. The key is to use Cauchy's inequality and distinct the cases of I = k and / = k + 1 
due to the weakest dissipative estimate of B. For I = k, we have 

2 

(2.63) 



v J V k u ■ V x V k B < e 
for I = k + 1, integrating by parts, we obtain 



V k+1 B 


2 


V k u 









L 2 



v I V k+1 u ■ V x V K+1 B = v I V x V K+1 u ■ V K+1 B < e 



7 k+l 



7k+l„ 



jk+l 



V k+1 B 


2 

+ c £ 


V k+2 u 




L 2 





L 2 



. (2.64) 
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Plugging the estimates ([23Tj) - (pTMj) and (I2371) - (l238|) from Step 1 into ^M\\ . by Cauchy's 
inequality, we then obtain 



d_ 

dt 



k+l 

£ 

l=k 

< e 



k+l 



[ V l u-V l E + \J2jv l E 
■ ■ i=k 

k+2 



V k+1 B 



2 

L 2 



Z=fc 



2 
2 

L 2 



(2.65) 



+ C k 



n, u 



'fc+2 



L 2 



+ 



V fc+1 B 



L 2 



This completes the dissipative estimate for E. 
Step 3: Dissipative estimate of B. 

Applying V fc to (|1.3P g and then taking the I? inner product with —V x V k B, we obtain 



V k d t E ■ V x V k B + 



V x V k B 



L 2 



< V 



V k u 



L 2 



V x V"B 



L 2 



+ V 



V k (f(n)u) 



L 2 



V x V k B 



(2.66) 



L 2 



Integrating by parts for both the t- and x-variables and using the equation (jl.3|) 4 . we have 



V k d t E ■ V x V k B = - — / V k E ■ V x V k B + I V x V k E ■ V k d t B 



dt 

= -. *L f V k E-V xV k B-v 
dtj 

From the estimates of I4 in Lemma 12.81 we have that 



V x V K E 



2 

L 2 



(2.67) 



V fc (/(n)n) r <C k \\u\\ HhnH2 V k (n,u) 



L 2 



(2.68) 



Plugging the estimates (|2.67p - (|2.68p into (|2.66p and by Cauchy's inequality, since divi? = 0, 
we then obtain 



d_ 

di 

< C 



V k E ■ V k V x B + A 



V k+l B 



L 2 



V k u 



L 2 



+ c 



V k+1 E 



L 2 



+ C k \\u\\ H k r\H 2 



V fe (n,u) 



(2.69) 



L 2 



This completes the dissipative estimate for B. 
Step 4-' Conclusion. 

Multiplying (|2.69p by a small enough but fixed constant rj and then adding it with (|2.65p so 
that the second term on the right-hand side of (|2.69p can be absorbed, then choosing e small 
enough so that the first term in (|2.65p can be absorbed; we obtain 



d_ 

dt 



/k+l 



l ^ J V l u ■ V l E - r, J V k E • V fc V xfijlA ^ || V l E 



\l=k 

k+2 

<cj2\\v l i 

i=k 



L 2 



+ 



V k+1 B 



L 2 



+ c k 



[n, u 



'k+2 



+ l|VB||i 2 £ V'(n, 



u 



J=k 



L 2 



+ 



L 2 



V k+1 B 



Adding the inequality above with (|2.59p . we get (|2.50p . 



L 2 I 
(2.70) 

□ 



2.3. Negative Sobolev estimates. In this subsection, we will derive the evolution of the 
negative Sobolev norms of (u,E,B). In order to estimate the nonlinear terms, we need to 
restrict ourselves to that s £ (0,3/2). We will establish the following lemma. 
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Lemma 2.10. For s G (0,1/2], we have 

d_ 
di 

and for s G (1/2,3/2), we have 



u,E,B)\\ 2 A _ s +X\\u\\%_ s < ([|(n,«)||^ + [|V5[|^J \\(u,E,B)\\ 6 _ s + \\E\\ 2 H2 ; (2.71) 



±\\( u ,E,B)\\%_ s +\\\u\\l_ s 

< (\\(n,u)f Hl + ||B||^ 1/2 \\VB\\% 2 ~ S \\u\\ L2 ) \\(u,E,B)\\ A _ s + 



(2.72) 



Proof. The A s (s > 0) energy estimate of (|1.3p 2 ()1.3j) 1 yield 

l£\\( u ,E,B)\\%_ s +u\\u\\%_ s 

= - j A~ s (u ■ Vu + unVn + u x B) ■ A~ s u + v J A~ s (f(n)u) ■ A' S E - J A~ s Vn • A~ s u 
< \\u ■ Vu + unVn + u x S||^_ s ||w||^- s + ||/(«HI/r- s ||-£||#- s + || Vn[|jj_ a ||tt||^_ a • 

(2.73) 

We now restrict the value of s in order to estimate the other terms on the right-hand side of 
(12T73L If s G (0,1/2], then 1/2 + s/3 < 1 and 3/s > 6. Then applying Lemma l2~H together 
with Holder's, Sobolev's and Young's inequalities, we obtain 

II" • Vn||^_ s < ||u • VuH^^. < ||«|| i3 /a ||Vn|| L2 

< ||Vn||^ 2 2+s ||V 2 n||^- s \\Vu\\v (2.74) 



< \\Vu\\ 2 H i + \\Vu\\ 2 L 2 . 



Similarly, we can bound 



||nVn||^_ s < ||Vn||^i + ||Vn|| 2 2 ; (2-75) 
\\uxB\\ A _ s < \\VB\\ 2 m + ||u|| 2 2 ; (2.76) 
\\f(n)u\\^ s <\\Vu\\ 2 Hl + \\n\\ 2 L2 . (2.77) 

Now if s G (1/2, 3/2), we shall estimate the right-hand side of (|2.73p in a different way. Since 
s G (1/2,3/2), we have that 1/2 + s/3 < 1 and 2 < 3/s < 6. Then applying Lemma 12.41 and 
using (different) Sobolev's inequality, we have 

||u • Vn||^_ s < \\u\\ l3/s ||Vn|| L 2 < ||u||^ 1/2 ||Vu||£ 2 2_s ||Vu|| i2 

< + WVufv ; ( 2 -78) 

\\nVn\\ A . s < ||n||^x + ||Vn|| 2 2 ; (2.79) 

lh x B\\h- s < WBW 6 ^ 2 \\VB\\% 2 - 3 \\u\\ L2 ; (2.80) 

\\f(n)u\\ k . s < \\u\\ 2 Hl + \\n\\ 2 L2 . (2.81) 

Note that we fail to estimate the remaining last term on the right-hand side of (|2.73|) as 
above. To overcome this obstacle, the key point is to make full use of the equation (|1.3j) 5 to 
rewrite n = n — /(re) + v~ 1 d\vE. This idea was also used in [181 121] • Indeed, using (jl.3p 5 . we 
have 

l|Vn[|a_. < ||A- s Vdiv£|| L2 + ||V(/(n) - n)\\ ti . s 
<ll^ + l|V(/(n)-n)||^ s . 

Here we have used the facts that s < 3/2 and f(n) — n = 0(n 2 ). Estimating the last term in 
(|2.82|) as before, and then collecting all the estimates we have derived, by Cauchy's inequality, 
we deduce lETTl) for s G (0, 1/2] and (I2T72D for s G (1/2, 3/2). □ 
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2.4. Negative Besov estimates. In this section, we will derive the evolution of the negative 
Besov norms of (u, E, B). The argument is similar to the previous subsection. 



Lemma 2.11. For s £ (0,1/2], we have 



!ll<«.s.B)ll5s 



+ A lid! 



< 



M||^ + ||VB||^ 



+ \\E\ 



H 2 i 



and for s £ (1/2,3/2], we have 



+ A IM 



< 



+ \\B\\" Ifi l, '\\VB\\fr \\u U2 



is-l/2 

II 2 



.3/2— s 

II 2 



(2.83) 



(2.84) 



z.oo 



+ ll£ll 



Proof. The Aj energy estimates of (|1.3|) 2 - (|1-3|) 1 yield, with multiplication of 2 2sj and then 
taking the supremum over j G Z, 



+ f kt 



< 



sup 2 



-2sj 



Aj (u • Vii + finVn + u x B) ■ Aju 



+ sup2- 2 ^' 
< || rz • Vn + jtinVn + u x B|| ^ 



i/ / Aj(f(n)u) ■ AjE - I AjVn ■ Aju 



u \\b- s + \\f( n ) u \ 



Br, 



\E\\ n-s + Vn Q- s 



(2.85) 

Then the proof is exactly same as the proof of Lemma [2.101 except that we should apply Lemma 
12.51 instead to estimate the B^^ norm. Note that we allow s = 3/2. □ 

3. Proof of theorems 

3.1. Proof of Theorem 11.11 In this subsection, we will prove the unique global solution to 
the system (|1.3p . and the key point is that we only assume the H 3 norm of initial data is small. 
Step 1. Global small £3 solution. 

We first close the energy estimates at the H 3 level by assuming a priori that \f £3 (t) < 5 is 
sufficiently small. Taking k = 0, 1 in (|2.20p of Lemma 12.81 and then summing up, we obtain 

3 3 



d_ 

lit 



L 2 



< V^3 + V^SV^S V^3 < SV 3 . 



1=0 1=0 
Taking k = 0, 1 in (|2.50p of Lemma 12.91 and then summing up, we obtain 



(3.1) 



JAY, /Vu-VVW^ fv l u-V l E- V J2 /VB-V'VxB] 

\i=o 1=0 1=0 J 

( 3 11 2 2 „ 2 2 M 

+ A E v ' n 2 +E v ^ 2 +E \\ vlB 



J=0 
3 



1=0 



1=1 



L 2 



(3.2) 



< 



1=0 



L 2 



+ 5 2 V 3 . 



Multiplying (|3.2|) by a sufficiently small but fixed factor e and then adding it with (|3.ip , since 5 
is small, we deduce that there exists an instant energy functional £3 equivalent to £3 such that 



-£3 + V3 < 0. 

at 



(3.3) 



Integrating the inequality above directly in time, we obtain (|1.10p . By a standard continuity 
argument, we then close the a priori estimates if we assume at initial time that £3(0) < 5q is 
sufficiently small. This concludes the unique global small £3 solution. 
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Step 2. Global £jy solution. 

We shall prove this by an induction on N > 3. By (jl.lOp . then (jl.lip is valid for N = 3. 
Assume (TOT]) holds for N- 1 (then now N > 4). Taking k = 0, . . . , N - 2 in (T2T20]) of Lemma 
and then summing up, we obtain 



d_ 

dl 



v 



^||V l (n ) u,£7,B) ^ 2 + A X]|| Viu 



v 



i=0 



2 

I? 



(3.4) 



Here we have used the fact that 3 < + 2 < A — 2 + 1 since A > 4. Note that it is important 
that we have put the two first factors in (|2.2U|) into the dissipation. 

Taking k = 0, . . . , N — 2 in (|2,5U[) of Lemma 12.91 and then summing up, we obtain 

'N-i . N-l r N-2 . \ 

/ V z n • VV'n + ^ / VV V l E ~vJ2 V/jE ' V x 



. 1=0 



1=0 



1=0 



N 



+ » E F 



d=0 

N 



9 Af-l 

+ v 

L 2 ^ 

/=0 



9 w-i 

+ v 

1=1 



V l B 



L 2 



(3.5) 



z=o 



<C^||v'u ^ + C N y/V N - iy /v^y/£^. 



Multiplying (|3.5p by a sufficiently small factor e and then adding it with (|3.4p . we deduce that 
there exists an instant energy functional £^ equivalent to £n such that, by Cauchy's inequality, 



This implies 



d 



— £n + — CnDn-iEn. 



(3.6) 
(3.7) 



(3.8) 



We then use the standard Gronwall lemma and the induction hypothesis to deduce that 

rt t 

£ N (t) + / V N {r)dT < C£ N (0)e CN foVN-i(r)dT 
Jo 

< C£ N {0)e CNPN - l{£N - l{0)) 

< C£ N (0)e CNPN -i {£N{0)) = P N (£ N (0)) ■ 

This concludes the global £n solution. The proof of Theorem 11.11 is completed. □ 

3.2. Proof of Theorem 11.21 In this subsection, we will prove the various time decay rates of 
the unique global solution to the system (jl.3p obtained in Theorem 11.11 Fix A > 5. We need 
to assume that £tv(0) < 5o = <5o(A) is small. Then Theorem 11.11 implies that there exists a 
unique global £n solution, and £/v(t) < Pn (£/v(0)) < 5o is small for all time t. Since now our 
5o is relative small with respect to A, we just ignore the A dependence of the constants in the 
energy estimates in the previous section. 
Step 1. Basic decay. 

For the convenience of presentations, we define a family of energy functionals and the corre- 
sponding dissipation rates with minimum derivative counts as 



ck+2 

c k 



k+2 

£ 

l=k 



\X7 l {n,u,E,B) 



L- 



and 



V 



k+2 



k+2 

EH 



k+l 



+E v lE 



l=k 



L 2 



+ 



V k+l B 



2 

L 2 



(3.9) 



(3.10) 
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By Lemma 12.81 we have that for k = 0, . . . , N — 2, 

k+2 „ „ k+2 „ 



|x;|v'(n,«,E )jB ) ; 2 +a^||v 

we have that for k = 0, . . . , N — 2, 



By Lemma 



l=k 
jk+2. 



2 

L 2 



(3.11) 



<Va DT i, + [|(T» > u)[| £ .||v^(n,ti) l2 v^(£,s) 



L-' 



d 

~dt 



/fe+i ~ fc+i » ~ 

(2 / vi,u • vv ' n + Yl / v ' n ■ vZ£; ~ ^ / 

\l=k J Z=fc J 



V k E ■ V fe V x 5 



/k+2 




fc+1 


2 






E 

\l=k 


V'n 


L 2 ^ II 


+ 

L 2 


V k+1 B 


1.) 



(3.12) 



fc+2 2 fc+2 

<E|| V ^ L2 + <5oE]| v W) 



2 

L 2 



Multiplying (|3.12p by a sufficiently small but fixed factor e and then adding it with (|3.1ip , since 
5o is small, we deduce that there exists an instant energy functional £ k+2 equivalent to £ k +2 
such that 



j/ k k +2 +V k k +2 <\\(n,u)\\ Lo 



V k+2 (n,u) ^ V k+2 (E,B) 



L 2 



(3.13) 



Note that we can not absorb the right-hand side of (|3.13p by the dissipation V^ 2 since it does 
not contain ||V fc+2 (-E, S)||^ 2 . We will distinct the arguments by the value of A;. If k = or 
k = 1, we bound \\V k+2 (E,B)\\ L2 by the energy. Then we have that for k = 0, 1, 



(It 



£ k k +2 + V k + 2 < Jv k+2 Jv k+2 ^£- 3 < A^ +2 , 



which implies 



_^_?fc+2 

dt k 



+ V k+2 < 0. 



(3.15) 



If k > 2, we have to bound \\V k+2 (E, B)\\ L2 in term of \\V k+1 (E, B)\\ L3 since yV k+2 can not 
control ||(n,u)|| Loo . The key point is to use the regularity interpolation method developed in 
[T7] . By Lemma 12.11 we have 

L 2 



[n,u 



V k+Z (n,u) ^ V K+A (E,B] 



<ll(n^)ll L2 2fe 



V k (n,u) 



V k+2 (n,u) V k+1 (E,B) 

L 2 



L 2 



\V a (E,B) 



where a is defined by 



3 \ 3 5 

k + 2 = (k + 1) x (1-7T7 +ax — => a = -k + l. 

2k I Ik 6 



(3.16) 



(3.17) 



Hence, for k > 2, if JV > \k + 1 2<k< |(iV - 1), then by (I3TTdT) . we deduce from (I3TT31 



that 



" pfc+2 , ^fe+2 < /F~~<r>fc+2 < /T~ ^fe+2 



(3.18) 



which allow us to arrive at that for any integer with < k < ^(N — 1) (note that N — 2 > 
|(iV - 1) > 2 since iV > 5), we have 



d ck+2 , -ptk+2 
dt k + k 



< 0. 



(3.19) 



The fact that T> k+2 is weaker than £ k i+2 prevents the exponential decay of the solution. 
In order to effectively derive the decay rate from (|3.19p . we still manage to bound the missing 
terms in the energy, that is, ||V & -B|L 2 and II V k+2 (E, B)\\ , 2 in terms of £ k+2 in (|3.19p . We again 
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use the regularity interpolation method, but now we need to also do the Sobolev interpolation 
between the negative and positive Sobolev norms. Assuming for the moment that we have 
proved (|1.14p or (|1.15p . Using Lemma I2UI we have that for s > and k + s > 0, 



V k B 



<^ IliJll fc + 1 + s 

L 2 - 11 



1 fc + s 

fc+1 



V k+L B 



k + l + a 



L 2 



fc + S 



fc+l + S 

L 2 



Similarly, using Lemma 12.7} we have that for s > and k + s > 0, 



(3.20) 



V k B 



< II Rll fc + 1 + s 

J-t 2 , oo 



fc + s 



L 2 



fc+l + S 

L 2 



(3.21) 



On the other hand, for k + 2 < N, we have 



< 



L 2 



V k+1 {E,B) 



N-k-l ||_JV 



L 2 



\L 2 



V fe+1 (-B,-B) 



N-k-2 



N ~ k ~ 1 (3.22) 

L 2 



Then we deduce from (|3.19p that 



dt k + r k J 



< 0, 



(3.23) 



where i? = max j^^, N \_ 2 j ■ Solving this inequality directly, we obtain in particular that 



-1/0 



£ fc fc+2 W < i £fe +2 (°) +M <C (l + t)- 1/,? = Co(l + t)- 



■ min{fc+s,A r -fc-2} 



(3.24) 



Notice that (j3.24p holds also for k + s = Qoxk + 2 = N. So, if we want to obtain the optimal 
decay rate of the whole solution for the spatial derivatives of order k, we only need to assume 
./V large enough (for fixed k and s) so that k + s < N — k — 2. Thus we should require that 



N > max { k + 2, -k + 1, 2k + 2 + s } = 2k + 2 + s. 



(3.25) 



This proves the optimal decay (j!.16|) . 

Finally, we turn back to prove (|1.14p and ()1.15p . First, we prove (|1.14p by using Lemma l2.10l 
However, we are not able to prove them for all s € [0,3/2) at this moment. We must distinct 
the arguments by the value of s. First, for s 6 (0, 1/2], integrating (|2.7ip in time, by (jl.lOp we 
obtain that for s G (0, 1/2], 



(u,E,B)(t)\\%_ s < \\(u ,E ,B )\\%_ 3 + [ V 3 (t)(1 + \\(u,E,B)(t)\\ A -.) dr 

Jo 



<C 1+ sup \\(u,E,B)(t) 



0<T<t 



\H~ S 



(3.26) 



By Cauchy's inequality, this together with (jl.lOp gives ()1.14p for s G [0, 1/2] and thus verifies 
(Ofi]) for s G [0,1/2]. Next, we let s G (1/2,1). Observing that we have (u ,E ,B ) G H~ 1/2 
since H- s n L 2 C H- s ' for any s' G [0, s], we then deduce from what we have proved for (|1.16p 
with s = 1/2 that the following decay result holds: 



fc+l/2 

V fc (n,n,E,5)(t) ^ 2 <C (l + t) ~ for A; = 0,1. 



(3.27) 
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(3.28) 



Here, since we have required N > 5 and now s = 1/2, so we can have taken k = 1 in fjl. 16j) . 
Thus by (EOTl) . (fLlOl) and Holder's inequality, we deduce from (ETT2D that for s G (1/2, 1), 

\\(u,E,B)(t)\\%_ a < \\(u ,E ,B Q )\\%_ s + f V s (t) (1 + \\(u,E,B)(t)\\ A -.) dr 

Jo 

+ /^|B(r)||^ 1/2 ||VB(T)|| 3 /^ s v^M||(«,£;,B)(T)|^_ s dr 

JO 

<C fl+fl+ f\l + T)- 2 ^- s / 2 Ur\ sup ||(«,^ 1 S)(r)y_. > ) 

<C fl+ sup ||(«,£;,B)(r)|| A _.y 

Here we have used the fact s G (1/2, 1) so that the time integral in (|3.28p is finite. This gives 
JLlif for s G (1/2, 1) and thus verifies (fl~T6|) for s G (1/2, 1). Now let s G [1,3/2). We choose 
so so that s — 1/2 < sq < 1. Hence, (uo, Eq, Bq) G H~ s °. We then deduce from what we have 
proved for f)1.16j) with s = sq that the following decay result holds: 



V k (n,u,E,B)(t) H2 <C {l + t) 



k + SQ 



for k = 0, 1. 



(3.29) 



Here, since we have required N > 5 and now s = so < 1, so we can have taken k = 1 in ([1,161) , 
Thus by (|3.29p and Holder's inequality, we deduce from (|2.72j) that for s € [1,3/2), similarly as 
in (I3T28D . 



(u,E,B)(t)\\l_ s <C (l+(l+ f\l + r)-^ +3 / 2 - s Ur) sup \\(u,E,B)(r)\\ H . 
V V Jo J o<r<t 

<C (l+ sup \\(u,E,B)(t)\\ 6 -.) . 



(3.30) 



Here we have used the fact s — sq < 1/2 so that the time integral in (|3.30p is finite. This gives 
(Oil) for s G [1,3/2) and thus verifies <fTTTTT|> for s G [1,3/2). Note that (fl~l~5j) can be proved 
similarly except that we use instead Lemma 12.111 
Step 2. Further decay. 

We first prove (|1.17p and (|1.18p . First, noticing that vf{n) = divE, by ()1.16p and Lemma 
221 if N > 2k + 4 + s, then 



V fc n(t) 



L- 



V*/(n)(0 



L 2 



v fc+1 ^) <a,(i+r^ 

.L 2 



(3.31) 



where we have used n = / _1 (/(n)). 

Next, applying V fc to (jl.3P n. (|1.3P o and then multiplying the resulting identities by X7 k u, X7 k E 
respectively, summing up and integrating over M 3 , we obtain 



1 d 

2dt 



V k (u,E) 



+ v 



V k u 



2 

L 2 



< 



J V fc (Vn + u-Vu + uriVn + u x B) ■ V k u + v J V k (V x B + f(n)u) ■ V k E 
V fc (u ■ Vu + fmVn + u x B) 



V fc+1 n 




V fc u 




L 2 







V k u 


L 2 





(3.32) 



+ 



V fc (VxB + f(n)u) 



L 2 



V k E 



L 2 



On the other hand, taking I = k in ()2.60p . we may have 

2 

L 2 



y V k d t u ■ V k E + ^ 

<( 



V k E 



V k+1 n 


+ 


V k u 


») 


V k E 


+ 




L 2 








L 2 



V fc (u ■ Vu + /mVn + u x B] 



x B) 




V k E 






L 2 




L 2 



(3.33) 
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Substituting (|2,6ip with I = k into (|3.33|) . we may then have 



4- I V k u-V k E~r 
dt 



V k E 



2 

L 2 



< c 


V k u 


2 / 

+ c( 


V k+1 n 


+ 


V k u 








L 2 V 




L 2 







+ 



V fc (n • Vn + unVn + n x B) 





V k E 




L 2 




L 2 



V k E 



+ 



L 2 



V k (V x B + /(n)u) 





V k u 


L 2 





L 2 

(3.34) 



Multiplying (|3.34p by a sufficiently small but fixed factor e and then adding it with ()3.32j) . 
since £ is small, we deduce that there exists equivalent to ||V fe (n, E 1 )^) L 2 such that, by 

Cauchy's inequality, ([2351 . (12391 . (l2Ti4l . (fl~16]l and (pOll . 



-J- fc (t) + J" fe (t) 
dt 



< 



< 



V k+1 n 



+ 



yfc+l^g 



+ 



V fc+1 (n,5) 



L 2 



L 2 



• " "if? 



V (n ■ Vn + unVn + n x 5) 



L 2 



+ 



V fc /(n)« 



L 2 



IV-BI 



L 2 



V k+l B 



L 2 



+ Hn,u) 



V k+1 (n,u) 



2 

L 2 



< C (l + t)-( k+1+s \ 



(3.35) 



where we required N > 2k + 4 + s. Applying the standard Gronwall lemma to (|3.35|) . we obtain 
W) < J- fc (0)e-* + C o f 6^-^(1 + T)^ k+1+s Ur <C (l + t)- (fc+1+s) . (3.36) 

JO 



This implies 



V fc (n,£)(t) <7^)<C7 (l + t)- 

L 2 



k + l+s 
2 



(3.37) 



We thus complete the proof of (|1.17|) . Notice that (jl. 18|) now follows by (|3.31|) with the improved 
decay rate of E in (|1.17|) . just requiring N > 2k + 6 + s. 

Now we prove (|1.19p . Assuming B^ = 0, then we can extract the following system from 
(|1.3|) 2 - (|1.3|) 3 . denoting ip = divn, 



dtn + ip = —u ■ Vn — /indivn, 

dtip + v4> — v 2 n = —An — div(n • Vn + unVn + n x B) + u 2 (f(n) — n). 



(3.38) 



Applying V to ()3.38p and then multiplying the resulting identities by v V n, V tp, respectively, 
summing up and integrating over 1R 3 , we obtain 



1 d 
2dt 



V k n 



+ 



v fc v 



+ v 



V fc ^ 



L 2 



-v 2 j V fc (n ■ Vn + undivu)V k n - J V k AnV k ^ 
J V k [div(n • Vn + unVn + n x B) - v 2 (f(n) - n)] V k tf). 



(3.39) 



Applying V k to (j3.38p 2 and then multiplying by — V fc n, as before integrating by parts over t 
and x variables and using the equation ()3.38p i . we may obtain 



d_ 

di 



V k i)V k n + v 2 



V k n 



2 

L 2 



L 2 



+ v / V fc nV fc V> + / V fc (n • Vn + /indivn)V fc V 



(3.40) 



J V fe [An + div(n • Vn + unVn + n x B) - v 2 (f(n) - n)] V k n. 
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Multiplying (|3.40p by a sufficiently small but fixed factor e and then adding it with (|3.39p . 

since e is small, we deduce that there exists Gk(t) equivalent to ||V fc (n, -0)|| L 2 such that, by 
Cauchy's inequality, 



j t G k {t) + g k {t) < 



V k+2 n 



V k+1 {u x B) 



L- 



+ 



V k+l (u-Vu) + V fe+1 (nVn) 



+ 



L 2 

V k (u-Vn) 



L 2 



+ 



V fc (ndivu) 



(3.41) 



L 2 



Notice that we have used the following estimates, by using the arguments in Lemma 12.2 



V fe (/(n)-n) 



L 2 



< c, lid 



ff3 



V k n 



L- 



<6n 



V k n 



2 

L 2 



(3.42) 



By Lemma [231 and Cauchy's inequality, we obtain 

2 



V fe+1 (-u x B) 



L 2 



u x V k+1 B + 



V k+ \u 



x B 



< 



u x V^ 1 ^ 



i2 
\L° 



L 2 



+ 



V k+ \u 



2 

L 2 

x B 



2 

L 2 



(3.43) 



V k+1 B 


2 2 

+ Vll roo 


V k B 


2 

+ 


V k+1 u 


2 




L 2 UL 




L 2 







l-BI 



The other nonlinear terms on the right-hand side of (|3.4ip can be estimated similarly. Hence, 
we deduce from (ET4"T|) that, by (fTTTUft — ffTTTST> . 



d_ 
dt 



Qk{t) + Q k {t) 



< 



V fc+2 n 


2 . 11 II 2 




2 

+ 




L 2 11 Hi 




L 2 



+ llVn| 



L°° 



+ \\(n,u) 



IL°° 



V fc+2 (n,n) i2 + ||V(n,«)|| 2 io 



V k B 



rk+li 



2 

+ 

L 2 


1 Rll 2 
\ n \\L°° 


V k+1 u 


(n,u) 


2 




L 2 





2 

L 2 



(3.44) 



< C ((1 + i)~( fe + 3 + s ) + (1 + t )-{k+V2+2s) + + ^-(fc+U/2+2.) 

< C (l + t)~ (fc+3+s) , 

where we required iV > 2/c + 6 + s. Applying the Gronwall lemma to (|3.44p again, we obtain 



Qk(t) < ^(O)e-* + Co f e^- T \l + r )-( fc+3+s ) dr < C (l + t)' 



(fc+3+s) 



This implies 



V*(n,^)(t) r2 < V / ^)<C (l + t)" 



fc+3+s 
2 



If required N > 2k + 10 + s, then by (|3.46|) . we have 



V k+2 n(t) r <C (l + i) 



fc+5+s 
2 



Having obtained such faster decay, we can then improve (|3.44p to be 
d Q k (t) + g k (t) < Co f (1 + t)-( fc+5+s ) + (1 + t)-(fc+r/a+2-) , \ < Co(1 + 1) 



Applying the Gronwall lemma again, we obtain 

<A/^y<c (i+r (fc/2+7/4+s) - 

L 2 



-(fc+7/2+2s) 



V*(ra,*/0(t) 



(3.45) 



(3.46) 



(3.47) 



(3.48) 



(3.49) 



We thus complete the proof of p,19|) . The proof of Theorem 11.21 is completed. 
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